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Abstract 

The C-parameter event-shape distribution for annihilation into hadrons is computed in the framework of 
SCET including input from fixed-order perturbation theory. We calculate all missing ingredients for achieving N^LL 
resummation accuracy in the cross section, which is then matched onto 0{a\) fixed-order results. Hadronization 
power corrections are incorporated as a convolution with a nonperturbative shape function. Wide-angle soft radiation 
effects introduce an (9(Aqcd) renormalon ambiguity in the cross section, which we cure by switching to the Rgap 
short-distance scheme. We also include hadron mass effects, but find their effect is rather small. Performing fits to 
the tail of the C-parameter distribution for many center of mass energies we find that the strong coupling constant is 
as(mz) = 0.1123 + 0.0015, with ;KVdof = 0.99. 
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1. Introduction 

The LEP collider, previously located at CERN, 
has delivered an enormous amount of highly accurate 
experimental data, which can be used to explore the the¬ 
ory of strong interactions in its high-energy regime. To 
study Quantum Chromodynamics (or QCD) at high en¬ 
ergies one needs to deal with jets: highly boosted and 
collimated bunches of particles that can be seen as the 
remnants of the underlying partons created at very short 
distances. One appealing strategy for describing jet dy¬ 
namics is through event shapes, infrared- and collinear- 
safe observables which are constructed from the energy 
and momenta of all the produced hadrons (in this sense 
event shapes are global quantities). They are designed 
to measure geometrical properties of the distribution of 
particles, and in particular they quantify how “jetty” the 
final state is. Additionally, being global observables, it 
is possible to compute high-order perturbative correc¬ 


tions, carry out log resummation to higher order, show 
factorization and exponentiation properties, and use fac¬ 
torization to control power corrections. 

One of the main uses of event-shape distributions is 
the determination of the strong coupling constant a^. 
The advantage of event shapes over other inclusive ob¬ 
servables is that they are essentially proportional to a^, 
rather than probing via corrections to a leading term 
(as is the case, for example, of DIS or the total hadronic 
cross section). Thus, event shapes are very sensitive to 
the strong coupling constant. On the other hand, event 
shapes are afflicted by nonperturbative power correc¬ 
tions and by large double Sudakov logarithms, which 
necessitate resummation. 

Here we study the C-parameter which can be written 

asEia: 

^ 3 S/jlAllRjIsin^ey 

2 (2,|p-|)^ • 
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It is interesting to compare C-parameter with thrust ||3l, 


T = 


-T = 


min 

Ft 



Pi\ \ 

2./ \Pj\ I ’ 


( 2 ) 


where n is referred to as the thrust axis. The three 
main differences between C and t are: a) C does not 
require a minimization procedure for its computation, 
whereas r does (namely hnding the thrust axis, event by 
event); b) C is dehned through a double sum, whereas 
T sums only particle by particle; c) the hxed-order 
prediction of C-parameter develops an integrable sin¬ 
gularity at Cshouider = 0.75, whereas thrust is always 
smooth. By shoulder we refer to the fact that the par- 
tonic cross section attains an integrable singularity Q at 
Cshouider, and only non-planar conhgurations contribute 
for C > Cshouider- 0 There are also a number of similari¬ 
ties between C and t, and perhaps the most remarkable 
one is that in the dijet limit (C, t « 1) and up to and in¬ 
cluding NLL resummation, both partonic cross sections 
are related in a simple way, which can be schematically 
expressed as follows: tnll = Cnll/6 0. Some other 
similarities will be highlighted later. 

Previous analyses of the thrust distribution using 
SCET at N^LL and analytic power corrections have 
found rather small (albeit precise) values of 
Two motivations for carrying out this new analysis are 
providing an additional determination of ctj and study¬ 
ing the universality of the leading power correction be¬ 
tween thrust and C-parameter. In this proceedings we 
summarize work done in Refs. nniiii]. 


2. Theoretical developments 

Until a few years ago, theoretical uncertainties were 
larger than the corresponding experimental ones and 
hadronic power corrections were not understood from 
ab-initio QCD considerations. The situation on the the¬ 
ory side has dramatically changed with the advent of 
Soft-Collinear Effective Theory (SCET) IIT2l IT3l [T4l 
[B). This effective held theory separates the relevant 
physics occurring at the various scales which play a 
role when jets are being produced: hard scale /lih, of 
the order of the center of mass energy Q (describes 
the production of partons at very short distances), jet 
scale Pj ~ Q VC/6 (describes the formation and evo¬ 
lution of jets at intermediate energies), and the soft 


*In Ref. (4) it is shown that soft gluon resummation at Cshouider 
makes up for a smooth distribution at LL order. 

^Other lower-order resummation event-shape analyses have also 
found small (although less accurate) values of a,, din- 


scale ps ~ QCI6 (describes wide angle soft radiation 
at longer distances). All three scales are widely sepa¬ 
rated for C « 1, and there is one function associated 
to each one of them: the hard coefficient Hq (the mod¬ 
ulus square of the QCD to SCET matching coefficient), 
which is common to all event-shape factorization theo¬ 
rems ; the Jet function Jj (built up with collinear Wilson 
lines), which is common for thrust, C-parameter na 
and Heavy Jet Mass (p) m ; and the Soft function Sc 
(dehned through soft Wilson lines), which in general de¬ 
pends on the specihc form of the event shape. Whereas 
the former two are perturbative (pH,pj » Aqcd), per¬ 
mitting the calculation of the hard and jet functions in 
powers of as, the soft function also has nonperturbative 
corrections that need to be accounted for (p^ ^qcd)- 
Renormalization evolution among the three scales sums 
up large logarithms to all orders in perturbation theory. 
It turns out that the anomalous dimensions for the soft 
function for C and t are identical, see cni. 

The soft function can be further factorized into a par- 
tonic soft function S^, calculable in perturbation the¬ 
ory, and a nonperturbative shape function Fc, which 
has to be obtained from hts to data[3 The treatment of 
hadronic power corrections greatly simplihes in the tail 
of the distribution, dehned by QCjb » Aqcd, where the 
shape function can be expanded in an OPE. The leading 
power correction is parametrized by Q^, the hrst mo¬ 
ment of the shape function. Interestingly, if one ignores 
hadron mass effects Eoiini, this matrix element can 
be related to the corresponding one in thrust in a trivial 
manner: 0[/2 = D.^li3n) = Qi 1221 . The main effect of 
this leading power correction is a shift of the cross sec¬ 
tion, dd-(C) —» dd-(C - Fl^lQ). When presenting the 
results of our hts, we will employ the power correction 
parameter Qi to ease comparison. 

The leading SCET factorization for the partonic 
C-parameter distribution can be written as EUIII]: 

1 do's Q r ' / QC s \ 

^ ^Ha(Q,p)jdsJAs,p)S,[— - -,p]. (3) 

It describes the most singular (and numerically dom¬ 
inating) partonic contributions in the dijet limit. The 
resummation of large logarithms is achieved by evolv¬ 
ing the functions Hq, J^, and from their respective 
natural scales pu, pj and ps, where logs are small, 
to a common scale p (which without loss of general¬ 
ity can be chosen to be, for instance, pj). In Eq. 0 
is also in the MS scheme, and suffers from an 


^Power corrections for the C-parameter distributions have been 
studied in other frameworks, see e.g. Refs. GDIIlllIS)- 
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(5(Aqcd) renormalon. We can switch to the renormalon- 
free Rgap scheme ^^2^ by performing subtractions on 
the partonic soft function (through an exponential of a 
derivative operator) and simultaneously allow the same 
terms to change from the MS scheme to the Rgap 
scheme. Our strong coupling will always be in the 
MS scheme. Adding these subtractions plus the renor¬ 
malization group evolution kernels gives; 

— ^ = ^HQ(Q,/iH)UH(Q,lJ^H,/i)[dsds'dk (4) 
o-Q dC 6 J 

X Jj{s,iij)U}{s - s',fd,fij)U^(k,ii,iis) 

a .. (QC - 371 A(R, Lis) s 

where 6(R,ids) is a series in powers of as(j-is) that can be 
computed directly from the thrust partonic soft function 
in Fourier space. For the renormalon to be properly can¬ 
celed by the subtractions, it is crucial that the exponen¬ 
tial and the partonic soft function are consistently ex¬ 
panded out to a given order in as{/is)- The subtractions 
introduce a scale R, which is close to the soft scale /ig 
and can be used to sum up large logs in the subtraction 
series through the finite shift parameter A{R,/is)- The 
dependence on R formally cancels between 6(R,ids) and 
A(R,/is) order by order in perturbation theory, but the R 
parameter is crucial to eliminate the Aqcd renormalon. 
Nonperturbative corrections are incorporated though a 
convolution with the shape function Fdp) whose first 
moment is The hadron level prediction for the dis¬ 
tribution is 



1 dcr 

o-Q dC 
dcr 
dC 





(5) 


and also includes the nonsingular contributions, 
d(Tns/dC, which in the dijet limit contains all terms 
which are kinematically suppressed by additional pow¬ 
ers of C. 

For our analysis we include perturbative corrections 
to the matrix elements Hq, and Sp to 0{al). For 
Hq they are known analytically, whereas for Jr and 
only the logarithmic terms at 0(al) are known (since 
the anomalous dimensions are known at three loops). 
These non-logarithmic terms are added as unknown co¬ 
efficients that are varied when we estimate the theory 
uncertainties. At 0{al) the jet function can be directly 
taken from Ref. ESll . The soft function needs to be 
computed to (9(aJ) ifTTl . which can be done analyti¬ 
cally at 0{as) and for the logarithmic corrections at 


For the non-logarithmic (9(a^) terms our eval¬ 
uation uses numerical output of the parton level MC 
EVENT-2 Il26ll27l . 

Through RGE evolution we achieve N^LL resumma¬ 
tion of the logarithmic terms. The anomalous dimen¬ 
sions required for solving the running equations can be 
taken directly from thrust. The only missing term is 
the four-loop cusp anomalous dimension, which is esti¬ 
mated using Fade approximants but is nevertheless var¬ 
ied in a wide range when estimating perturbative un¬ 
certainties. Its effect is in any case negligibly small. 
The required components for a given resummation or¬ 
der are specified in Table We introduce a primed 
counting, which is defined as the regular (unprimed) 
one, but with the matrix elements being included to one 
order higher. Eor consistency, the renormalon subtrac¬ 
tion series are including to one order higher as well. The 
primed counting achieves a better description of data 
and allows the correct summation of logs at the level of 
the distribution (for an extended discussion of this the 
reader is referred to 128]). 

We include nonsingular terms at the same order as 
the functions Hq, Jr, S^. These can be obtained by 
subtracting the fixed-order singular cross section as de¬ 
scribed by the SCET factorization theorem from the 
full QCD partonic distribution. The latter can be com¬ 
puted analytically at 0(as), and determined numerically 
at Oia^) and 0(al) from the parton-level MC programs 
EVENT2 and EERAD3 ll29l l30ll . respectively. Eor the 
0(a^) and 0(a^J nonsingular contributions our numeri¬ 
cal procedure entails uncertainties which are accounted 
in the estimate of the theoretical uncertainty. 

It is sometimes customary to write the most singular 
terms of an event-shape cumulant cross section in the 
following exponentiated form: 


2 (C) = - 

o-Jo 


0-Jo dC' 



X exp 


' oo /+! 

V/=l j=l 


as(Q) \ 
Itt 


lnO| 


(6) 


Erom the result for the factorization theorem in Eq. 0 
one can determine the Gij and coefficients as shown 
in Table 1^ see ifTOl . 


Resummation Order 

Calculable G;/s and B/s 

EL 

C /,,+1 

NEE' 

G/,/ and Bi 

N^LL' 

Gi,/-I and B 2 

N^LL' 

G /,/_2 and B 3 
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cusp 

non-cusp 

matching 

P 

ns 

6 

LL 

1 

- 

tree 

1 

- 

- 

NLL 

2 

1 

tree 

2 

- 

- 

N^LL 

3 

2 

1 

3 

1 

1 

N^LL 

4P 

3 

2 

4 

2 

2 

NLL' 

2 

1 

1 

2 

1 

1 

N^LL' 

3 

2 

2 

3 

2 

2 

N^LL' 

4P 

3 

3 

4 

3 

3 


Table 1: Loop corrections for primed and unprimed orders. For the 
anomalous dimensions of A(R,ps) one uses the same orders as for 
other non-cusp anomalous dimensions. The superscript “p” indicates 
that a Fade approximation is being used. 

3. Setting the Renormalization Scales 

The C-parameter can be divided into three distinct re¬ 
gions, in which the renormalization scales must satisfy 
different constraints 

1) nonperturbative: C < Stt Aqcd 

I^H ~ Q, ^J^J ~ V^QCd2> rs Aqcd , 

2 ) resummation; Stt Aqcd C < 0.75 (7) 

[C QC 

~ Q, ^J^J ~ Q A/ -7, -^ » Aqcd , 

Vo o 

3) fixed-order: C > 0.75 

- Mj = - R ~ Aqcd ■ 

These three regions are sometimes referred to as the 
peak, tail and far-tail regions, respectively. In order to 
satisfy these requirements we need to use renormaliza¬ 
tion scales that depend on the value of C, called profile 
functions. The constraints in Eq. Q do not fully specify 
the profile functions, but this ambiguity cancels order- 
by-order in perturbation theory. This allows variations 
of the profiles to be used to estimate perturbative uncer¬ 
tainties. The specific form of the profile functions and 
the variation of their parameters are given in Ref. IfTTII . 
and illustrated in Fig. [T] 

In Fig. we show our C-parameter cross section for 
Q - niz, together with experimental data. This hgure is 
produced with our best theoretical prediction and uses 
our central values for afiniz) and Qi presented in Sec.|^ 
The center blue line corresponds to the prediction for 
our central profiles, whereas the light blue band shows 
the perturbative uncertainty. 

4. Fit results 

Fitting for Qi together with as(inz) accounts for 
hadronization effects in a model-independent way. In 



Figure 1: Profile functions for the renormalization scales /in, 
and fis(C), when using the default profile parameters (center thick 
line) and when varying them (light band). Fully canonical profiles are 
shown in gray. 


o-dC 



Figure 2: Theoretical prediction for the C-parameter distribution at 
N^LL’ order for Q = mz, using the best fit values for as(mz) and Hi. 
The blue band corresponds to the theory uncertainty as described in 
the text. Experimental data for various experiments are also shown. 


order to determine afiniz) and Qi in the same fit, one 
needs to perform a global analysis that includes data at 
many center of mass energies Q. For each Q the differ¬ 
ential distribution has a noticeable degeneracy between 
the two fit parameters, and the use of data from the dif¬ 
ferent Q values breaks the degeneracy. Hence LEP and 
SLAC data are employed together with data from lower 
energy experiments such as TRISTAN and PETRA. For 
our analysis we use all available experimental data with 
energies between 35 GeV and 207 GeV in the tail re¬ 
gion. To estimate theoretical uncertainties we perform 
500 fits at each order in the resummation, NLL', N^LL', 
and N^LL', with theory parameters randomly chosen for 
each fit. These parameters specify; the profile func¬ 
tions, unknown perturbative coefficients, or statistical 
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I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 

I N'^LL' + Geo's) + Q.[(Rp) 
■ N-LU + 0(0'^) + ni(R,/j) 



0.1000 0.1025 0.1050 0.1075 0.1100 0.1125 0.1150 


asimz) 

(a) 



ctsimz) 

(b) 


Figure 3: The left panel (a) shows the distribution of best fit points in the a^imz) - 2Q,i plane for fits performed with our best theoretical predictions: 
resummation of large logs and power corrections defined in the Rgap scheme with renormalon subtractions. The dashed lines corresponds to an 
ellipse fit to the contour of the best-fit points to determine the theoretical uncertainty. The total (experimental + theoretical) 39% CL standard 
eiTor ellipses are displayed (solid lines), which coiTespond to l-cr (68% CL) for either one-dimensional projection. The big points represent the 
central values in the random scan for cusimz) and 2 Hi. The right panel (b) shows the distribution of best fit points in the cts{mz) -^^/dof plane, 
corresponding to the points given in panel (a). 


order cTjOwz) (with Qi) q'^(otz) (with Qi(/;a,/^a)) 
NLL' 0.1071(60)(05) 0.1059(62)(05) 

NNLL' 0.1102(32)(09) 0.1100(33)(06) 

N^LL' 0.1117(16)(06) 0.1123(14)(06) 

Table 2: Best fit values for asimz) at various orders with theory uncer¬ 
tainties from the theory scan (first value in brackets), and experimen¬ 
tal and hadronic eiTor added in quadrature (second value in brackets). 
Our final resulf at N^LL' is shown in bold face. 

uncertainties on the numerical determination of the non¬ 
singular contributions. The result of these many fits are 
shown graphically as dots in Fig. We show two pro¬ 
jections; as vs 2Q.I in panel (a), and vs ^min/dof 
in panel (b). As the resummation order increases the 
perturbative uncertainty decreases as expected, and the 
also decreases signihcantly. The correspond¬ 
ing numerical results and uncertainties are shown in Ta- 

bleHim. 

In Fig. I^we show determinations of Usiniz) from hts 
to the C-parameter distribution with different levels of 
theoretical sophistication. From left to right they are: 
fixed order with 0{a\) (large logs not yet summed up), 
N^LL' resummation (no Qi in the fit), with power cor¬ 
rections included (not yet removing the renormalon). 


including Rgap scheme (not yet accounting for hadron 
masses), and with hadron mass effects. All error bars 
are perturbative, so the error bars of the hrst two deter¬ 
minations to the left of the vertical dashed line do not 
account for the neglect of power corrections. Including 
the nonperturbative power corrections by fitting has 
the greatest effect on asimz). Hadron mass effects give 
negligible contributions within current uncertainties. 


5. Conclusions 

We have presented an accurate determination of 
asimz) from hts to the C-parameter distribution. The 
key points to our precise theoretical prediction are: a) 
higher order resummation accuracy (N^LL), achieved 
through the SCET factorization theorem, b) inclusion of 
Oial) matrix elements and hxed-order kinematic power 
corrections, c) held-theoretical treatment of nonpertur¬ 
bative power corrections, and d) switching to a short- 
distance Rgap scheme, in which the sensitivity to in¬ 
frared physics is reduced. We have not discussed hadron 
mass effects, as their effect is quite small. A thorough 
discussion can be found in IfTOll . 
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a^imz) from global C-parameter tail fits 
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_ 
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1 

Below error bars & ± -> perturbative error 

- 

Cfs(mz) 


oca's) fixed-order i 

All errors: ttsfmz) = 0.1123 ± 0.0015 

- 



' 0.1317 ±0.0052 1 
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1 

+ R-scheme ^ . 



1 
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" 
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1 
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" 
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Figure 4: Impact of the different components of our theoretical setup on the determination of a^(mz). 


Our final results from the global analysis reads im 

asimz) = 0.1123 + (0.0002)exp (8) 

+ (0.0007)hadr ± (0.0014)pert , 

= 0.1123 + (0.0015)tot 

We conclude by presenting a comparison of our 
C-parameter fit with the determinations of ctj and Qi 
from our previous thrust analysis © of the thrust dis¬ 
tribution, see Fig. The figure shows that the de¬ 
termination of the strong coupling constant for both 
event shapes is compatible. Moreover there is uni¬ 
versality in the result for the leading power correction 
= Q^/(3;7r) = f^j/2 determined in both analyses. 
The two independent determinations agree within their 
1 -cr uncertainties, where the precision of the extrac¬ 
tions is greater than that achieved in the past. Note 
that without including the respective prefactors (37r and 
2) (shown in green) the values of Qj = 0.329 GeV 
and = 1.98 GeV numerically differ by a; 4.5-cr, so 
the agreement nicely demonstrates the consistency of 
our theoretical predictions. A detailed comparison with 
other Usimz) determinations can be found in Ref. ca. 
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